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1. Introduction 

The Density Matrix Renormalization Group (DMRG) as a major technique for 
one-dimensional quantum systems provides a variational method in the space of 
matrix-product states . The key idea of DMRG comprises the repeated basis 
truncation by density matrix projection in an iteratively enlarged system. This idea 
turned out to be extraordinary successful beyond the original purpose of computing 
the low-energy spectrum of quantum chains with short-range interactions 

Finite-temperature properties of quantum chains can be obtained by map- 
ping the one-dimensional quantu rn system to a two-dimensional classical lattice 
via Trotter-Suzuki decomposition EEMinilll, The Transfer Matrix DMRG method 
(TMRG) ^'^ yields the finite temperature properties of infinite-sized quantum 

chains by applying a DMRG algorithm to the transfer matrix of such a correspond- 
ing two-dimensional lattice in order to find the highest contributing eigenvalues. 

In this paper, we propose a different method: The Quantum Corner Transfer 
Matrix DMRG ( QCTMRG) for the calculation of thermodynamic properties of finite 
one-dimensional quantum systems is a variant of Nishino's Corner Transfer Matrix 
DMRG (CTMRG) for classical systems ^i^. Here, both calculating the partition 
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sum and obtaining the reduced density matrix is just one step and does not involve 
finding eigenvectors of large transfer matrices. Thus, CTMRG performs drastically 
better than the TMRG method for classical systems. Can we benefit from the ideas 
of the CTMRG in the quantum case? 

In contrast to the classical models treated with CTMRG so far, two substantial 
differences arise in dealing with a Trotter decomposition: 

• A significant anisotropy due to the existence of a well-distinguished real 
space and a Trotter direction. 

• Calculation of the trace in order to obtain the partition function demands 
periodic boundary conditions of the two-dimensional plane. 

Both aspects will be considered our QCTMRG approach. In this paper, we present 
the QCTMRG algorithm and its essential features and close with first results and 
a discussion of the method's performance. 

2. CTMRG for quantum systems 
2.1. Trotter decomposition 

We consider one-dimensional quantum chains of length L where the Hamiltonian 
H = X^j^i ^j.j+i only includes nearest-neighbor interactions. Generically, the par- 
tition function of such a quantum chain does not factorize because the neighboring 

interaction terms do not commute. This difficulty has been overcome by the Trotter- 
. . Is) . 

Suzuki decomposition ^ which maps the partition function of the quantum chain 
onto a partition function of a two-dimension al c lassical chequerboard model. A 
variant was introduced by Sirker and Kliimper where the partition function 

Z = tre-^^= lim tr \T^e~^"ne-'"]^'^^ (1) 

where e = P/M, is decomposed into a trace over a product of imaginary-time 
propagators. Here, Tr and Tl are left and right shift operators. Note that because 
of translational invariance, [Tr^l,^^] = 0. The decomposition with shift operators 
overcomes disadvantages of the common cheq uerboard decomposition because the 
spatial periodicity of the lattice is one site^^^. Yet, the QCTMRG algorithm as de- 
scribed in the subsequent section can be adapted to the chequerboard decomposition 
with only minor changes. 

As a benefit from the factorization, we are able to insert identities of the form 

^|si...Si)(si...Si| = 1 (2) 

pictured as slices of discrete temperature or imaginary time, where the chain state 
|si . . .sl) is the tensor product of local quantum states \sj) at site j. Thus, we 
achieve a classical two dimensional model, spanned by the real space or chain di- 
rection and the auxiliary introduced imaginary time or so called Trotter direction. 
The local transfer matrix (see figure [1]) 
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Fig. 1. Graphical representation of 4-spin transfer- matrix t. The bended curves illustrate the 
symmetries of the object. 



{siSi+i\T\siSi+i) := {siSi+i\e ^'^'■''+^\siSi+i) (3) 

is a fourth order tensor of dimension n x n x n x n where n is the number of states 
per site. With this, we find 

L L 

(si . . . sije^^-^lsi . . . sl) = ^ J]^ (^iSi+i|'^|'SiMi+i) = X! H (si-iMiKl^i-iSi)- (4) 

j = l.-.L j=l...L 

With the action of the right and left shift operators 

L 

{s'l... sl|7r,l|si ...sl)= Y]_^'>i±i,S'' (^) 

i 

we get the partition function 

L M/2 



of a two-dimensional classical lattice. For finite M the corrections to the approxi- 
mated partitions functions and free energies have been shown to be no larger than 
©(e^) li^. Graphically (see figure [2]), the decomposition is represented as rows of 
alternately rotated local transfer matrices r. 

Thermal expectation values of local operators can be calculated by a modified 
partition function. We recall the statistical definition of the thermodynamical ex- 
pectation value 

m = ^^ with Z{0,)=trO,e-^" (7) 

of an operator O ai a site i. In terms of a two-dimensional classical model, the non- 
normalized expectation value Z(Oi) is formed by replacing one standard transfer- 
matrix within the product ([6]) by the modified local transfer-matrix f := 0^6"*^''. 
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real space direction L 



Fig. 2. Graphical representation of the Trotter decomposition 

2.2. Constituting tiles 

We consider the partition function of a Trotter decomposition Starting with 
an initial lattice small enough to be treated analytically, our aim is to iteratively 
expand this system in both directions to large enough sizes. Periodic boundary con- 
ditions correspond to the trace in the calculation of the partition function. They 
are physically essential in Trotter direction. We are, however, free to choose open 
boundary conditions in real space direction. This choice significantly reduces com- 
putational effort because the tensor-dimensionality of the corner tiles is now three 
rather than four in the case of real-space periodic boundary conditions. Thus, free 
sites at both real space edges will be integrated out, whereas the edge states re- 
garding Trotter direction are left as degrees of freedom in order to permit periodical 
closing. 

At an arbitrary renormalization step, we consider a lattice of fixed size in both 
directions divided into four parts (see figure[3]). Assuming free edges in Trotter direc- 
tion, we no longer deal with two-dimensional matrices as in conventional CTMRG, 
but with three-dimensional tensors. So the pieces of the system are third and fourth 
order tensors and, thus, shall rather be called tiles than-somewhat misleadingly- 
matrices in this context. Here, a composition of left and right corner tiles Cl,Cr 
makes up the whole system. 

Now we introduce some additional tiles which are required for enlarging the 
system in the renormalization step (figure |4]) . The unit cell of the two-dimensional 
Trotter decomposition forms the smallest basic tile W with four free edges. Arising 
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Fig. 3. Trotter decomposition segmented into four corner tiles. The periodical boundary condi- 
tions are illustrated by the double-headed arrow. The corner tiles Cl,Cii have been depicted by 
special symbols reflecting the underlying symmetry. 
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Fig. 4. Graphical representation of left and right corner tiles CljCji, row-to-row tiles Pl,Pji, 
column-to-column tile Py and basic tile W. An additional tile X is used for the measurement of 
expectation values. 



from the spatial symmetries, three tiles play the role of the row-to-row transfer 
matrix of CTMRG. Here, we have for both left and right side the row-to-row tiles 
Pl and Pr and a single column-to-column tile Py, which are tensors of third resp. 
fourth ordeiQ. 



''Depending of the symmetries of the spin chain and the Trotter decomposition, Cl and Cr as 
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Together with these tiles, we have to keep account of the different bases asso- 
ciated with their edges. The basic tile W demands a vertical basis {|cry)} at the 
left and right edge and a horizontal basis {|T/f^} at the upper and lower edge, 
both consisting of one- or two-site states depending on the size of the unit cell of 
the Trotter decomposition. These bases will be unaffected by changes within the 
renormalization procedure. The edge bases of the corner tiles, however, are itera- 
tively renormalized during the CTMRG-iterations to keep a fixed maximum size. 
The maximum size of the vertical-edge basis {jwy)} of the corner tiles is my. The 
horizontal edges of the corner tiles are different for the left/right corner tile. Their 
bases {|tol)}, {j'Ttfl)} have a maximum size of rriL resp. ttir. Tiles Pl, Pr and Pv 
include bases which stay unaltered as well as renormalized bases corresponding to 
their different edges. 

2.3. Renormalization algorithm 

With the concept of the constituting tiles, we sketch the outline of the QCTMRG 
algorithm: 

(1) Construction of initial tiles 

In the special case of the Sirker-like decomposition the basic tile W (figure[S]) 
as well as the column-to-column tile Py are built by addition of two mutually 
rotated transfer matrices r. The left row-to-row tile Pl and the left corner tile 




S 

Fig. 5. Graphical representation of the basic tile W in the Sirker-decomposition. A composition 
of two transfer matrices can be depicted as a tile with four arms at the edges. The shape of the 
arms symbolizes the symmetries of the tile and the basis of the underlying tensor operator. 

Cl arise by summing out the left free sites of W. For the right row-to-row tile 

well as and Pr can be represented by a single tile in special cases. Here, we assume the general 
case with no vertical reflection symmetry. 
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Pr and the right corner tile Ca the right free sites of W are bent into the 
horizontal edges in order to aim a summation with the neighboring site on the 
lower/upper tile when composing tiles. In figure [6l the construction of all initial 
tiles is depicted. These initial tiles demand the construction of the initial basis. 
Starting from the one-site basis {jc)}, the vertical basis of corner and column- 
to-column tile i^ ~ {["'v)} = {I"") ® I"')}' while the horizontal basis of 
left corner and row-to-row tile is {|™l)} = — {|'^)}' ^hc horizontal 
basis of right corner and row-to-row tile is {jm^fl)} — {jc) ® |(t)} in the first 
iteration. 




Fig. 6. Graphical representation of the initial tiles. 

(2) Calculation of expectation values 

All tiles are combined to form a periodically closed two-dimensional lattice in 
order to determine the partition function of the system. Summing out the states 
on the inner edges, we obtain the desired partition sum Z of the system. The 
partition sum Z{0) of the modified system with a certain operator O situated 
in the middle of the lattice (see figure [7]) can be realized similarly. Thus, the 
thermodynamical expectation value (O) = Z{0)/Z can be computed. 

(3) Enlargement of system tiles 

In the next step of the renormalization process, the system is expanded by 
enlargement of the corner tiles in applying a row-to-row, a column-to-column 
tile and the basic tile. Correspondingly, the row-to-row and column-to-column 
tiles have to be expanded by addition of the basic tile. These enlargement steps, 
which technically correspond to matrix-matrix-multiplications, are illustrated in 
figurelHl Enlargement of the tiles implies enlargement of the bases, which is done 
by simple tensor products { = { |"^v)® lev)} , { I^^l)} = { jwi)® |cr_f/)} 
and { I m^;) } = { | ^r) (X" | ) } . Note that the sizes of the bases grow by a factor 
of 4 to 16 depending on the underlying spin system, which correspondingly 
increases the size of the tiles. 

(4) Construction of reduced density matrices 



'■(55 denotes the usual tensor product. 
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(a) 










Fig. 7. (a) Calculation of the partition function Z and (b) of the partition function of the system 
containing a certain operator Z{0). 
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Fig. 8. Enlargement of system tiles, (a) The enlarged corner tile Cj^ is composed of the corner tile 
Cl, the column-to-column tile Py, the row-to-row tile and the basic tile W. (b) The enlarged 
the column-to-column tile Py is composed of the column-to-column tile Py and the basic tile W. 



The crucial step in the DMRG-like renormalization procedure is the construc- 
tion of the reduced density matrix. In the Trotter decomposed lattice, three 
density matrices pv, Pl and pn are required in order to renormalize the three 
different types of bases appearing in the system (see figure [5]). The concept of 
the reduced density matrix is to compose the corner matrices and to sum out 
all but one edge. The pictu re of the reduced density matrices then corresponds 
to cuts of the system ^iSIIHl 
(5) Truncation of basis 

The idea behind renormalization group procedures is to iteratively integrate 
out insignificant degrees of freedom. In the context of DMRG-type algorithms. 



April 11, 2008 10:57 WSPC/INSTRUCTION FILE Bartel-schad2 



Quantum C orner- Transfer Matrix DM RG 9 




Fig. 9. Construction of the (a) vertical and (b) left reduced density matrix. Illustration as "cuts 
of the system" . 

measuring the contribution of states for calculating the partition function is 
carried out by diagonalization of the reduced density matrix. Those eigenstates 
with large eigenvalues will dominate because the partition function is nothing 
but the trace over the reduced density matrix. Conservation of the bases' sizes 
demands truncation of the eigenstates with lowest weight. So, we establish a 
projection onto the m states with largest eigenvalues as renormalization pre- 
scription. Since a projection is obtained for each of the three bases, all edges of 
the system tiles are now reduced to a fixed size keeping only the most relevant 
states for calculation of the partition function. 
(6) Iteration 

Go to step ([2]) until desired system size is reached. 

2.4. Normalization of growing tiles 

While our interest lies in the calculation of (local) expectation values of a certain 
quantum mechanical system, we have to deal with the partition functions of an it- 
eratively increased classical system in the QCTMRG-algorithm. Thus, the partition 
function is a rapidly growing entity leading to several huge matrix entries in the 
tiles' numerical representation in each renormalization step. So, the program runs 
the risk of exceeding the numerical capacity of the variables of the system. To avoid 
this problem, a constant prefactor is extracted in each renormalization step. 

2.5. Implementation 

Our implementation of the QCTMRG method follows the algorithm presented in 
section 12.31 The ratio e = /3 /M is kept at a fixed value during the renormalization 
procedure as in TMRG. This provides, in contrast to TMRG, not only a new (lower) 
temperature but also a different (larger) chain length after each renormalization 
step. Measurements of expectation values can be made in the center of odd chains 
or at the two sites in the center of even chains. Our implementation makes use of 
good quantum numbers which drastically reduce computational effort. 
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Fig. 10. Free energy per site vs. temperature for tlie (classical) antiferromagnetic Ising chain of 
varying chain length. The free energy (circles) from QCTMRG with m = 32 and e = 0.05 agrees 
to high precision to the exact value (full line). Note that the chain length is L = 10/T. The free 
energy of the infinite chain is plotted (dashed line) for comparison. 



3. Results 

As a first test, we conipute the energy and the free energy of the Ising model 
The antiferromagnetic Ising chain (J > 0) with the classical Hamiltonian 

i/ = J^a,a,+i (8) 

3 

is exactly solvable by a transfer matrix method The classical Ising spin a can 
take the values +1,-1. 

We have computed the thermodynamics of the Ising chain using the QCTMRG 
algorithm, keeping m = 32 states within the renormalization procedure. The inverse 
factor of temperature T and number of imaginary time steps M has been chosen as 
e — (TM)^^ — 0.05. We expect an excellent agreement with the analytical results 
because the Trotter decomposition becomes exact for the classical Ising model. In 
figures [TOl and [TT] the numerical results are plotted and, indeed, both data show the 
expected agreement. 

We face a different situation in considering the antiferromagnetic spin-1/2 
Heisenberg chain 

H = JJ2S,-S,+,. (9) 

i 

The model is exactly solvable but quantum fluctuations are no longer suppressed. 
Thus, the Heisenberg chain can serve as a non-trivial trial system for the QCTMRG 
algorithm. 
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Fig. 11. Expected local energy in the center of the chain vs. temperature for the antiferromagnetic 
Ising chain of varying chain length. The local energy expectation value in the center of the chain 
(circles and squares denote odd and even chain lengths) from QCTMRG with m = 32 and e = 0.05 
agrees to high precision to the exact local energy (full line). 
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Fig. 12. Free energy density vs. temperature for the Heisenberg chain of varying chain length. The 
free energy density is plotted for different QCTMRG calculations (stolid lines) with a preserved 
number of states m = 50 (circles), 150 (squares), 400 (diamonds). Note that the chain length is 
L = 10/T. The free energy density (data from TMRG, m = 200) of the infinite chain is plotted 
(dotted line) for comparison. 



We calculated the free energy and the expectation value of the energy operator 
in the middle of the chain, see Figs. [121 [HI The calculations have been done with 
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fixed ratio e = 0.05 while the preserved number of states m from the renormalization 
was varied from m = 50 to m = 400. 

In figure [T2l the free energy density is plotted against temperature. Note that 
the chain length is related to temperature by L = lO/T. With increasing m the free 
energy tends to converge. Interestingly, the convergence is faster for lower values 
of T and correspondingly larger system sizes. We account for that point later in 
Section 14.21 For a comparison, we added the well-converged data for an infinite 
chain calculated by conventional TMRG. We see the deviation of the QCTMRG 
data from this curve more pronounced at higher temperatures which are related to 
smaller system sizes. Thus, we can interpret the deviation as finite size effect. 
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Fig. 13. Free energy density vs. inverse number of states for the Heisenberg chain of various 
lengths at various temperatures. The free energy densities of chains with T = 0.25J and L = 
40 (diamonds), T = O.IJ and L = 100 (squares), and T = 0.04J and L = 250 (circles) have 
been calculated by QCTMRG (open symbols) with a varying number of states kept during the 
renormalization. The dashed lines are linear fits. The filled symbols denote the free energy of an 
infinite system calculated by TMRG. 



For three points of fixed temperature and chain length the free energy density 
was extrapolated, see figure [TS] We find / growing with 1/m for each point. This 
results from the fact that the partition sum is underestimated for smaller number 
m of states kept. This leads to a monotonically increasing free energy for growing 
1/m. As already noticed earlier, we find better convergence for smaller temperature 
and larger system sizes. 

The energy expectation value in the center of the chain has been calculated for 
various system sizes and temperatures in figure [T4l Again we have e = 0.05 which 
leads to a fixed relation LT = 10 for the QCTMRG data. Odd and even chain 
lengths are included in this calculation, in contrast to the free energy data which 
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Fig. 14. Energy per site vs. temperature for the Heisenberg cliain of varying chain length. The 
energy expectation value at the chain center is plotted for different QCTMRG calculations with 
a preserved number of states m = 50 (circles), 150 (squares), 400 (diamonds). The upper points 
belong to systems with an odd number of sites, the lower points belong to systems with an even 
number of sites. Note that the chain length is L = 10/T. The exact energy density of the infinite 
chain [22] is plotted as a full line for comparison. 



have only been given for even chains. For odd chains, the expectation value was 
taken from a plaquette at the exact center of the decomposition. For even chains, 
we considered one of the two central plaquettes. Even chain lengths appear as the 
natural choice in the renormalization procedure. Chains with an odd number of 
sites have a well-defined central site which might give better results. 

We find the energy expectation value to lie for even systems below and for odd 
systems above the infinite-system value. The low-energy spectrum of the antiferro- 
magnetic Heisenberg chain involves spin- 1/2 spinous which can be identified with 
quantum domain walls. The ground state is a total spin singlet on chains with an 
even number of spins. In odd chains there is no spin singlet ground state and so 
always a spinon "excitation" is present. In the domain wall picture, there is always 
a kink present. For this reason, we expect the chains with an even number of sites 
to possess a lower local energy than chains with an odd number of sites. However, 
we expect that the expectation values for even and odd chains converge to the same 
limit in an infinite chain. This description agrees well with our observed behaviour. 

In either case, we have a strong dependence of the expectation value on the 
number of preserved states m in the renormalization group. For both even and odd 
number of sites, we can still distinguish the data curves up to low temperatures 
even for high values of m. Yet, the convergence for even system sizes is faster than 
for odd sizes. 
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Fig. 15. Energy per site vs. inverse number of states for the Heisenberg chain of various lengths at 
various temperatures. The energy expectation values at the center of even chains (lower symbols) 
with T = 0.25J and L = 40 (diamonds), T = O.IJ and L = 100 (squares), and T = 0.04J and 
L = 250 (circles) and of odd chains (upper symbols) with T = 0.256 J and L = 39 (diamonds), 
T = O.IOIJ and L = 99 (squares), and T = 0.0402J and L = 249 (circles) have been calculated. 
We denote the exact energy density of the infinite chain by full symbols. The dashed lines are 
linear fits. 



To get a more quantitative picture we plotted the convergence of several points 
of fixed temperature and system size, see figure [13 Like in the free energy data, the 
high-T and larger sized systems show a better convergence. We included linear fits 
for a better understanding. The expectation values for even/odd sized system at 
finite temperatures seem to differ even in an 1/m extrapolation. The infinite-chain 
expectation value lies well-between those boundaries. For larger system sizes and 
lower temperatures we obtain a better convergence in 1/m. 

The results from this section shall serve as an illustrating background for a 
discussion of the scope of application and the limits of the QCTMRG technique. 



4. Discussion of the QCTMRG method 

A crucial attribute for the QCTMRG algorithm is the fixed relation between system 
size and temperature during the renormalization loops. This condition sets the focus 
of the method to finite size systems, in contrast to the TMRG. Yet, a significant 
advantage in computational resources would make the QCTMRG an interesting 
choice for low-temperature studies of very large systems where finite-size effects 
play a marginal role. 
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4.1. Running time and storage use 

For a view on the algorithm's running time and storage use, we studied the scaling 
behaviour with the number of states kept m of the different steps in the algorithm. 
We find the algorithm to take asymptotically 0(771^) elementary floating point op- 
erations (FLOP) and a storage use of order 0{m^) floating point numbers for large 
m. However, most of the time is spent in matrix-matrix-multiplications operations 
where a naive approach would give a running time of 0{N^) FLOP for multipli- 
cation of two N X N matrices. In contrast to this, the fastest algorithm currently 
known has an asymptotic running-time of 0{N'^ '^'^^) FLOP^^Sl, So we expect the 
running-time to scale asymptotically with well-below fourth order in m, in a clever 
implementation. 

In case of the TMRG algorithm, the by-far most time consuming part is finding 
the largest eigenvalues of the transfer-matrix. We, thus, expect the asymptotic scal- 
ing behaviour to be dominated by the employed Arnoldi method which performs 
with between 0(1) and 0{N) x the cost of a certain matrix-vector product when 
N X N is the size of the matrix. In our case, the necessary matrix-vector product 
has a cost of 0{m^) FLOP and the total matrix has a basis of size N oc m^. We 
thus expect the total running-time to vary between ©(to^) worst case and 0{m^) 
best case. Further reduction might be achieved when clever algorithms for matrix 
multiplications are implemented. The storage scales with the size of the system- 
block-transfer-matrix (O(m^) numbers). 

From the theoretical point of view, the QCTMRG and the TMRG algorithms 
have similar running-times in the large-m limit. The TMRG algorithm will be favor- 
able in those cases where the Arnoldi method reaches a fast convergence. For some 
ill-posed problems, however, the QCTMRG algorithm might have an advantage. 
The TMRG algorithm is the clear winner when storage usage is a sensible quantity. 

In our implementation of both methods, we observed a roughly similar running 
time of both algorithms for different number of preserved states m. The storage use 
was, indeed, higher in the QCTMRG algorithm. 

4.2. Fundamental challenges 

Yet, we still face some fundamental challenges of the QCTMRG algorithm in ad- 
dition to the advantages of TMRG concerning the use of system resources. As 
already mentioned above, the quantum character of the system enforces periodical 
boundary conditions in Trotter direction which makes the two-dimensional corner 
transfer-matrix from Nishino's CTMRG a three-dimensional tensor in the case of 
QCTMRG. This is the reason behind an increase in running-time and storage which 
is crucial and partly destroys the advantages from the development of CTMRG over 
TMRG. 

Another more subtle point involves the spectra of the reduced density matrices 
from the renormalization procedure. Consider the limit T 00 which is close to 
the starting point of TMRG and QCTMRG where (3 = 1/T ^ eM is small. Here, 
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the local transfer-matrix of the Trotter decomposition reduces to 

{s'is'i+i\T\siSi+i) := (siSi+i|e"^'*''*+i|siSi+i) = dsi,siSsi+i,si+i (10) 

which means that the initial spin configuration will not be changed by the transfer- 
matrix. 

Now, we consider a Trotter decomposition for infinite temperature or vanishing e. 
The chosen graphical representation depicts the Kronecker symbols as lines passing 
through the transfer-matrix plaquettes. If the system is periodically closed in Trotter 
direction and has open boundaries in space direction, we can imagine the paths as 
non- interacting lines around a cylinder. This depicts the trace from the partition 
sum. 

Building a reduced density matrix introduces a vertical or horizontal cut into the 
system. Consider the case of a cut in Trotter direction which corresponds to the cut 
in TMRG. One "spin path" has been cut through while all the others stay intact, i.e. 
they still are summed out in the reduced density matrix. Only the intersected thread 
determines the eigenspectrum of the reduced density matrix. With a straightforward 
calculation we find one or two degenerate eigenvalues depending on whether we 
intersect an odd or even number of sites. All other eigenvalues remain zero in this 
case. This turns the renormalization step, which is a truncation of basis states, to be 
highly effective. In a system with small e we still find a small number of dominating 
eigenvalues which lead a DMRG renormalization to success. 

The opposite situation is faced when the cut is made in space direction. If the 
underlying spin chain has 2L spins, the density matrix will "cut" as much as L 
paths. All L spins act separately and, thus, all spin configurations have the same 
contribution to the partition function. Merely the remaining degrees of freedom will 
be summed out. Consequently, we are left with a reduced density matrix which has 

degenerate eigenstates when S is the spin size. No effective truncation can be 
found. The situation in cases with high, but finite temperature is certainly less ill- 
posed. Though, still we expect a slow decay of the spectrum of the reduced transfer 
matrix. 

This scenario explains why the calculated QCTMRG data shows a most sig- 
nificant deviation from the expected values at higher temperatures. Starting with 
a small system size the algorithm can still handle the system with high precision. 
For larger system sizes, a truncation has to be done during the renormalization 
step which cannot be optimal since the reduced density matrix will still have a 
flat spectrum. At even larger system sizes and lower temperatures the spectrum of 
horizontal density matrix will become acceptably well-behaved again. 

5. Conclusions and Outlook 

We developed a new method for finite-temperature studies of one-dimensional quan- 
tum systems based on the CTMRG of Nishino. The free energy densities and thermal 
energy expectation values at the chain centers have been successfully calculated for 
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the classical Ising chain and the antiferromagnetic spin- 1/2 Heisenberg chain by 
the Quantum Corner- Transfer Matrix DMRG. Reliable results were given for finite 
temperatures and system sizes. Yet, the algorithm faces two diflaculties: 

• Periodic boundary conditions reduce the efiiciency. 

• The reduced density matrix in space direction has a slowly decaying eigenspec- 
trum. 

If finite-temperature data in the thermodynamic limit are aimed at, the quantum 
TMRG method is certainly still the method of choice. At least one of the mentioned 
problems should be solved to make the QCTMRG technique an attractive option. 
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